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A  thin  layer,  Navier-Stokes  flow  solver  was  used  to  predict  the  aerodynamics  about 
interfering  bodies.  Surface  pressure  distributions  compared  very  well  with  avail¬ 
able  experimental  data  for  freestream  Mach  numbers  between  0.60  and  1.20  at  an 
interbody  separation  distance  of  0.8  diameters.  Five  separate  interference  regions 
were  identified  which  contribute  to  the  total  interference  force,  which  was  an 
attractive  force  for  all  cases  investigated.  In  general,  the  level  of  interference 
increases  with  increasing  free  stream  Mach  number,  and  decreases  with  increasing 
separation  distance. 
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ABSTRACT 

A  thin-layer,  Navier-Stokes  flow  solver  was  used 
to  predict  the  aerodynamics  about  interfering  bod¬ 
ies.  Surface  pressure  distributions  compared  very  well 
with  available  experimental  data  for  freestream  Mach 
numbers  between  0.60  and  1.20  at  an  interbody  sep¬ 
aration  distance  of  0.8  diameters.  Five  separate  in¬ 
terference  regions  were  identified  which  contribute  to 
the  total  interference  force,  which  was  an  attractive 
force  for  all  cases  investigated.  In  general,  the  level 
of  interference  increases  with  increasing  free  stream 
Mach  number,  and  decreases  with  increasing  separa- 
tion  distance. 

INTRODUCTION 

The  aerodynamic  properties  and  the  performance 
qualities  of  fighter  aircraft  are  influenced  by  the 
aerodynamic  forces  generated  by  mutually  interfering 
stores  in  external  carriage.  In  addition,  a  knowledge 
of  these  forces  is  necessary  for  determining  the  tra¬ 
jectories  of  these  bodies,  relative  to  each  other  and 
relative  to  the  aircraft,  as  they  are  released. 

Experimental  and  computational  investigations  of 
the  aerodynamics  of  mutually  interfering  bodies  have 
been  performed  for  single,  double  and  triple  combi¬ 
nations,  Refs  [1,2].  In  these  investigations,  Euler  so¬ 
lutions  were  obtained  and  compared  with  the  experi¬ 
mental  measurements. 

In  the  present  work,  viscous  flow  solutions  are  ob¬ 
tained  to  investigate  the  aerodynamics  of  mutually 
interfering  bodies.  A  main  element  of  the  study  is 
the  investigation  of  the  forces  imposed  on  the  bod¬ 
ies  when  they  are  in  close  proximity.  Previous  studies 
were  performed  for  a  fixed  distance  between  body  cen¬ 
terlines.  Each  of  the  studies  considered  a  range  of  sub¬ 
sonic  and  supersonic  free-stream  Mach  numbers.  The 
study  in  Ref  [2]  showed  that  the  finned  bodies  tend  to 
be  pulled  together  at  subsonic  free-stream  conditions 
and  to  be  repulsed  at  supersonic  free-stream  condi¬ 
tions.  However,  the  repulsive  forces  were  relatively 
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weak  in  comparison  to  the  attractive  forces.  The  im¬ 
portance  of  viscous  effects  increases  in  determining 
the  interactive  forces  between  the  bodies  as  the  dis¬ 
tance  decreases  between  them.  The  use  of  the  viscous 
solver  allows  us  to  investigate  the  interactive  forces 
acting  on  the  bodies  at  different  conditions  for  which 
the  use  of  the  inviscid  solver  would  not  be  valid. 

GRID  GENERATION 

The  body  configuration  is  a  generic  store  with  a 
cylindrical  centerbody,  and  a  tangent  ogive  forebody 
and  afterbody.  The  afterbody  is  truncated  to  mount 
to  a  sting.  The  radius  of  the  cylindrical  centerbody 
is  one  inch  and  the  total  length  of  the  store  is  11.882 
inches  from  the  nose  to  the  truncation  point  at  the 
sting.  The  tangent  ogive  forebody  and  afterbody  are 
given  by  circular  arcs  with  radius  of  6.056  inches.  The 
computational  grid  for  these  configurations  was  gen¬ 
erated  by  the  Program  EAGLE-Numerical  Grid  Gen¬ 
eration  system,  Ref  [3].  The  grid  generated  for  the 
two  and  three  store  configurations  is  an  elliptic,  five 
block,  O-C-H  grid  system.  The  grid  i$  a  single  ogive- 
cylinder-ogive  with  sting  and  utilizes  two  reflection 
planes  to  simulate  both  the  two  store  and  three  store 
conditions. 

A  block  structure  was  chosen  which  could  be  eas¬ 
ily  modified  to  change  the  separation  distance  be¬ 
tween  the  stores.  The  block  structure  is  basically 
the  same  for  the  two  and  three  store  configurations. 
A  cross-sectional  cut  of  the  grid,  at  a  constant  / 
plane,  is  shown  for  the  two  store  configuration  in 
Figure  1(a),  and  the  three  store  configuration  in  Fig¬ 
ure  1(b).  Blocks  1  and  3  are  embedded  C  grids  which 
surround  the  store  and  sting.  The  C  configuration  was 
chosen  to  maintain  orthogonality  of  grid  lines  off  the 
nose  of  the  ogive  forebody.  For  these  two  blocks,  the 
/  lines  of  the  computational  grid  run  longitudinally, 
along  the  body  surface,  the  J  lines  are  normal  to  the 
surface,  and  the  K  lines  are  circumferential  about 
the  body.  Blocks  2  and  4  are  immediately  in  front  of 
blocks  1  and  3.  The  /  lines  in  blocks  2  and  4  run 
longitudinally,  along  the  stagnation  line  out  the  nose 
of  the  body,  J  lines  are  normal  to  the  stagnation  line, 
and  K  lines  are  circumferential  about  the  stagnation 
line.  Transition  is  made  from  the  constant  1  planes  of 
’blocks  1  and  3  to  the  constant  J  planes  of  blocks  2  and 
4.  This  is  more  easily  seen  in  Figure  1(c),  which  shows 
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a  constant  K  plane  cut  through  the  grids.  Block  5 
transitions  the  inner  blocks  1  and  2,  along  the  reflec¬ 
tion  planes,  to  the  farfleld  boundaries.  The  I  lines 
are  longitudinal,  J  lines  are  normal,  and  K  lines  are 
circumferential.  Blocks  3  and  4  are  constrained  by  the 
reflection  planes  and,  therefore,  require  no  transition 
to  farfioM  boundaries. 

The  reflection  planes  were  aligned  along  the  planes 
of  symmetry  for  each  of  the  configurations.  The  in¬ 
cluded  angle  between  the  two  reflection  planes  of  the 
two  store  configuration  is  90  degrees,  while  that  of  the 
three  store  configuration  is  60  degrees.  For  purposes 
of  analysis,  the  angle  0  is  defined  as  positive  clock¬ 
wise,  with  0  =  0  starting  at  the  inboard,  constrained 
flow  location,  to  0  =  x  at  the  outboard,  unconstrained 
flow  location.  This  is  shown  in  Figure  1(a)  for  the  two 
store  case  and  Figure  1(b)  for  the  three  store  case. 

The  store  has  97  points  distributed  in  the  stream- 
wise  direction  with  concentration  of  the  surface  points 
at  the  expected  areas  of  aerodynamic  interest,  the 
ogive  nose,  the  shoulder  transition  to  the  cylindrical 
midsection,  and  the  ogive  boattail.  Additionally, the 
32  points  in  the  normal  direction  (J)  were  concen¬ 
trated  close  to  the  body  surface  and  smoothly  transi¬ 
tioned  to  the  outer  boundaries.  The  total  dimensions 
for  blocks  1  and  3  are  97  x  32  x  16  and  97  x  32  x  6, 
for  a  total  of  21  circumferential  points.  Dimensions 
of  the  remaining  blocks  were  20  x  20  x  16  for  block  2, 
20  x  20  x  6  for  block  4,  and  116  x  30  x  16  for  block  5. 

FLOW  SOLVER 

The  solutions  were  obtained  using  a  thin-layer, 
Navier-Stokes  solver  based  upon  that  of  Ref  [4].  The 
solver  uses  a  flux-difference  split  scheme  based  on 
Roe’s  approximate  Riemann  solver.  The  governing 
equations  are  written  in  strong  conservation  law  form 
to  capture  shock  waves  in  the  flow  field.  A  finite  vol¬ 
ume  formulation  is  used  to  discretize  the  flow  govern¬ 
ing  equations.  The  solver  uses  an  implicit,  two-pass 
scheme.  Local  time-stepping  is  used  to  accelerate  the 
convergence  of  the  solution.  Characteristic  variable 
boundary  conditions  are  used  at  the  farfield  bound¬ 
aries  and  no-slip  boundary  conditions  are  used  at  the 
body  surfaces.  A  multi-block,  body  conforming  com¬ 
putation  mesh  is  used  in  the  computations.  Turbu¬ 
lence  is  modeled  by  using  the  Baldwin- Lomax  model. 

RESULTS 

Computations  were  performed  for  the  flow  about 
interfering  bodies  at  dilferent  free  stream  Mach  num¬ 
bers,  Mao  i  and  separation  distances  D,  where  D 
is  the  distance  between  the  bodies  at  their  closest 
points.  A  Reynold’s  number  of  1.2  x  106  was  cho¬ 
sen,  so  that  the  quality  of  computational  results  ob¬ 
tained  could  be  measured  against  the  available  exper¬ 
imental  data  of  Ref  [5].  Extensive  comparisons  were 


made  for  the  varying  values  of  M »  and  circumferen¬ 
tial  body  locations  at  D  =  1.6  inches.  Figures  2(a), 
(b),  and  (c)  show  a  comparison  of  pressure  coefficient 
along  the  surface  for  the  two  body  configuration  at 
D  =  1.6  inches,  JlfM  =  0.95.  The  pressure  distribution 
for  the  most  constrained  flow,  the  inboard  location 
0  —  0,  is  shown  in  Figure  2(a).  Computational  pre¬ 
dictions  agree  very  well  with  the  experimental  results 
through  both  the  expansions  at  the  nose  and  boat- 
tail,  and  quite  accurately  predicts  the  shock  location 
and  strength  of  compression  at  mid-body.  A  slight  de¬ 
viation  in  the  prediction  of  the  minimum  Cp  is  noted. 
Pressure  distributions  at  0  =  j  and  0  =  x  are  shown 
in  Figures  2(b)  and  (c),  respectively.  In  general,  the 
computational  predictions  agree  well  with  a  slight  un¬ 
derprediction  of  the  shock  strength  and  recompres¬ 
sion.  A  random  selection  of  experimental  comparisons 
for  the  two  b  jdy  case  at  =  0.6,  0.85, 1.05,  and  1.2 
are  shown  in  Figures  3(a),  (b),  (c),  and  (d).  Again, 
the  predicted  Cp  values  agree  quite  well  with  exper¬ 
iment  for  these  cases.  It  should  be  noted,  however, 
that  for  the  M =  0.85  case  the  shock  location  is 
predicted  slightly  forward  of  the  experimental  result, 
and  the  strength  of  compression  is  underpredicted. 
To  some  extent,  this  deviation  can  be  attributed  to 
insufficient  mesh  concentration  in  this  area. 

Interference-free  flows,  for  the  case  of  a  single 
body,  are  axisymmetric.  The  introduction  of  a  neigh¬ 
boring  body  results  in  flow  conditions  which  are  differ¬ 
ent  in  the  unconstrained  flow  region  neighboring  the 
0  =  x  plane  and  the  constrained  flow  region  neigh¬ 
boring  the  0  =  0  plane.  This  destroys  the  flow  ax- 
isymmetry.  Therefore,  the  deviation  of  the  flow  from 
symmetry  about  the  xy  plane  and  the  deviation  of 
the  flow  from  axisymmetry  are  measures  of  flow  in¬ 
terference.  Examples  of  parameters  which  may  be 
used  as  measures  for  flow  interference  are  <rj,  <Ji  and 
A.  The  parameter  <j\  is  based  on  comparing  the  flow 
properties  along  the  body  at  the  minimum  interfer¬ 
ence  station,  0  —  ic,  and  the  maximum  interference 
station,  0  =  0.  It  is  given  by 

/  |ACp  (x,0)|  dx 

ai  = - 71 - 

where  A Cp  (x,  0) ,  A Cp  (x,  0)  =  Cp  (x,  x)  -  Cp  (x,  0) , 
is  the  difference  between  pressure  coefficient  distribu¬ 
tions  at  the  inboard  station,  0  =  0,  and  the  outboard 
station,  0  =  x.  The  integral  is  taken  over  the  body 
length  and  d  is  the  body  diameter.  The  ^  factor  ap¬ 
pearing  in  the  definition  of  a\  is  simply  a  scaling  factor 
to  allow  comparisons  between  <j\  and  a2,  which  is  de¬ 
fined  below.  Although  the  parameter  rrl  takes  into  ac¬ 
count  the  flow  properties  at  only  two  angular  stations, 
0  =  0,  and  0  =  x,  it  is  a  indication  of  interference  over 
the  total  surface,  since  the  flow  along  all  angular  po¬ 
sitions  is  dependent  on  the  conditions  at  0  =  0  and 
0  =  x.  The  parameter  crj  is  obtained  by  integrating 
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the  absolute  value  of  A Cp  (x,9)  over  the  surface  area 
projected  onto  the  xy  plane  for  —  ^  <  9  <  j,  where, 
ACP  (x,6)  =  Cv  (x,  x  —  9)  —  Cp  (x,0) .  Therefore, 

_  /  \ACp  ( x,6 )  et  •  raj  dS 

«(i)2 

where  ez  and  n  are  the  unit  vector  in  the  z  direction 
and  the  unit  outward  normal  to  the  body  surface,  re¬ 
spectively,  and  the  integral  is  taken  over  the  body 
surface.  The  parameter  A  is  based  on  measuring  the 
flow  deviation  from  axisymmetry.  It  is  given  by 

,  /  I  A/ 1  dS 

S 

where  5  is  the  body  surface  area,  and  A /  =  fi  — 
foy  where  /  is  some  flow  property.  Subscript  i  refers 
to  the  interference  flow,  while  subscript  o  refers  to  a 
reference  flow  distribution  which  may,  for  example,  be 
chosen  to  be  that  of  an  interference-free  flow  (single 
body),  or  some  average  distribution  obtained  from  the 
interference  flow  such  as 

, ,  .  _  CjiWydo 

lx 

The  values  of  the  different  parameters  suggested 
above  are  not  expected  to  be  equal,  however,  the  use 
of  any  one  of  them  allows  a  measure  of  interference  to 
be  obtained  for  comparing  the  flow  at  different  con¬ 
ditions.  The  parameters  a\  and  <r3  are  used  in  the 
present  study.  The  parameter  A,  however,  is  just  given 
as  an  additional  example  for  possible  interference  pa¬ 
rameters. 

The  interference  force  coefficient,  Ca,  acting  on  a 
body  due  to  a  neighboring  body  is  obtained  by  in¬ 
tegrating  ACp{x,$)  over  the  surface  area  projected 
onto  the  xy  plane  for  —  j  <  9  <  Positive  A Cp 
values  contribute  to  an  attractive  interference  force, 
while  negative  A Cp  values  contribute  to  a  repulsive 
interference  force.  Comparisons  were  made  along  the 
body  between  pressure  distributions  at  9  =  0  and 
0  =  x.  This  resulted  in  identifying  distinct  regions 
along  the  body  where  attractive  or  repulsive  interfer¬ 
ence  forces  act.  It  is  possible  to  identify  these  re¬ 
gions  on  the  basis  of  comparisons  at  the  0  =  0  and 
9  =  x  angular  positions,  since  pressure  integrals  at 
these  positions  are  the  most  highly  weighted  among 
all  angular  positions  when  computing  aerodynamic  in¬ 
terference  forces  and  since  the  flow  in  the  9  direction 
is  generally  continuous  and  gradual  variations  in  flow 
properties  take  place  as  9  varies  between  the  stations 
9  =  0  and  9  =  x.  The  total  effect,  whether  repul¬ 
sive  or  attractive,  depends  on  the  integral  effects  of 
all  regions. 

As  the  flow  moves  along  the  body,  it  passes 
through  regions  of  expansion,  compression,  expan¬ 
sion,  and  compression,  respeatively.  For  a  given  sepa¬ 
ration  distance,  D,  the  extent  of  each  of  these  regions 


along  the  body  depends  on  the  free-stream  Mach  num¬ 
ber,  Moo  •  The  upstream  and  downstream  boundary 
positions  for  each  of  these  regions  are  functions  of  the 
angle  9}  due  to  the  different  flow  conditions  existing  in 
the  unconstrained  flow  region  neighboring  the  9  =  x 
plane  and  the  constrained  flow  region  neighboring  the 
9  =  0  plane.  These  different  flow  conditions  also  cause 
the  strength  of  the  expansions  and  compressions  to 
depend  on  9. 

Figure  4(a)  compares  the  pressure  distribution 
along  the  inboard  body  stations  at  9  =  0  (curve  i)  and 
the  outboard  body  stations  at  9  =  x  (curve  o)  for  the 
two  body  case  at  Mx  ~  0.95  and  D  =  1.6  inches.  The 
nose  expansion  region,  shown  in  the  figure  depends  on 
Mqo  •  At  relatively  low  values,  it  ends  upstream 
of  the  front  shoulder.  The  downstream  boundary  of 
this  region  moves  toward  the  shoulder  as  the  free- 
stream  Mach  number  increases.  For  all  values 
considered  here,  this  region  does  not  extend  beyond 
the  shoulder  in  the  unconstrained  outboard  position, 
9  =  x.  However,  in  the  constrained,  inboard  posi¬ 
tion,  9  =  0,  it  continues  to  extend  further  beyond  the 
shoulder  as  M^  continues  to  increase.  The  x  position 
of  the  downstream  boundary  of  the  expansion  region 
is  generally  a  decreasing  function  of  9.  The  pressure 
in  this  region  is  higher  in  the  inboard  side  than  the 
outboard  side.  Due  to  the  relatively  early  termina¬ 
tion  of  this  region  in  the  outboard  side,  a  cross  over 
between  curves  i  and  o  occurs  at  point  B.  The  re¬ 
gion  extending  between  points  A  and  B  is  the  front 
repulsive,  FR,  region.  Downstream  of  point  B,  the 
flow  continues  to  expand  in  the  inboard  side  and  this 
expansion  is  finally  terminated  by  a  relatively  strong 
shock.  In  the  outboard  side,  the  flow  goes  through 
a  compression  region  terminated  by  a  relatively  weak 
shock.  The  shock  surface  terminating  the  compression 
region  moves  downstream  as  M<x,  increases.  However, 
it  moves  at  a  faster  rate  in  the  outboard  position  than 
at  the  inboard  position.  Therefore,  a  second  cross  over 
point  occurs  at  point  C,  creating  a  center  front  attrac¬ 
tive,  CFA,  region  between  points  B  and  C.  As  the  flow 
continues  to  compress  beyond  point  C,  then  expand 
and  compress,  a  center  repulsive,  CR,  region  between 
points  CandD,  a  center  back  attractive,  CBA,  region 
between  points  D  and  E,  and  a  back  repulsive,  BR, 
region  between  points  E  and  Fare  formed.  In  general, 
large  interference  forces  associated  with  the  individ¬ 
ual  regions  indicate  large  flow  interference,  however, 
the  combined  effects  of  the  different  regions  is  not  nec¬ 
essarily  proportionate  to  the  effects  of  the  individual 
regions. 

The  extents  and  effects  of  the  five  regions  identi¬ 
fied  above  are  strongly  dependent  on  M^,  .  Under  cer¬ 
tain  conditions,  some  of  these  regions  may  be  absent. 
Figure  4(b)  compares  the  inboard  pressure  distribu¬ 
tions  (curve  t)  and  the  outboard  pressure  distribution 
(curve  o)  for  the  two  body  case  at  MK  =  0.60  and 
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D  =  1.6  inches.  At  this  low  Moo  value,  the  flow  is 
nearly  axisymmetric  and  deviations  between  curves  t 
and  o  are  small.  Also,  the  flow  is  shockless.  The  FR 
region  is  apparent  in  this  figure,  however,  the  inter¬ 
ference  effects  are  not  sufficiently  strong  to  cause  the 
front  outboard  compression  region  to  extend  down¬ 
stream  of  the  front  inboard  compression  region.  As  a 
result,  the  CR  region  is  not  formed,  and  the  CFA  and 
CBA  regions  combine  to  from  a  single  region.  Fig¬ 
ure  4(c)  compares  the  pressure  distributions  at  the 
inboard  and  outboard  stations  for  Moo  —  1-2.  Com- 
paiing  this  figure  to  Figure  4(a),  it  is  apparent  that 
regions  CFA  and  CR  are  further  expanded  as  Moo  in¬ 
creases,  while  the  BR  region  is  completely  eliminated. 

The  pressure  comparisons  corresponding  to  the 
Moo  values  for  Figures  4(a),  (b),  and  (c)  are  shown 
in  Figures  5(a),  (b),  and  (c),  for  D  —  0.8  inches.  The 
force  regions  are  apparent  in  these  figures.  The  re¬ 
duced  separation  distance  affects  the  pressure  distri¬ 
bution  both  at  the  inboard  and  outboard  stations, 
however,  the  strongest  effects  are  apparent  at  the  in¬ 
board  stations,  where  higher  pressure  gradients  are 
observed. 

Figure  6(a)  compares  the  pressure  distributions  of 
the  two  body  case  at  Moo  =  0.95  for  the  three  separa- 
tion  distances  D  =  0.8,  1.2,  and  1.6  inches  at  C  =  0. 
A  similar  comparison  is  made  in  Figure  6(b)  at  6  =  rr. 
It  is  clear  from  Figure  6(a)  that  a  reduced  separation 
distance  causes  the  flow  to  reach  lower  pressure  lev¬ 
els  in  the  nose  expansion  region,  which  is  terminated 
with  a  stronger  shock  than  the  corresponding  shock 
at  larger  D  values,  and  which  is  positioned  upstream 
of  those  in  the  larger  Q  cases. 

Figure  7  compares  a\  and  oj  for  the  1.6  inch  sep¬ 
aration  distance,  while  a  similar  comparison  is  pre¬ 
sented  in  Figure  8  for  the  0.8  inch  separation  dis¬ 
tance.  It  is  clear  that  both  parameters  predict  the 
same  trends  in  interference  variations  with  Moo  ■  For 
the  large  separation  distance,  D  =  1.6  inches,  it  is 
observed  that  the  level  of  interference  increases  with 
Moo  ,  for  the  range  of  values  considered  here.  How¬ 
ever,  the  rate  of  increase  in  interference  levels  with 
Moo  declines  at  the  larger  Moo  values.  At  the  smaller 
separation  distance,  D  =  0.8  inches,  the  level  of  in¬ 
terference  increases  with  Moo  i  until  the  flow  becomes 
supersonic.  Beyond  that  point  the  interference  levels 
off  and  further  increases  in  Moo  cause  it  to  decline 
slightly. 

Figure  9  shows  the  variation  of  the  two  body  in¬ 
terference  force  with  free  stream  Mach  number,  for 
the  two  separation  distances  considered  here.  For 
the  Moo  range  considered,  the  interference  force  is 
an  attractive  force.  For  a  given  Moo  value,  the  at¬ 
tractive  force  increases  with  reduced  distance.  For 
low  Moo  values,  the  magnitude  of  the  force  increases 
with  Moo  until  it  reaches  its  maximum  value,  at  a  free 
stream  Mach  number  of  approximately  0.8.  Beyond 


that  point  the  attractive  force  is  reduced  with  in¬ 
creased  Moo  ,  except  for  the  local  maximum  at  about 
a  Moo  value  of  1.05.  Computations  for  the  separation 
distance  of  1.2  inches  were  limited  to  the  Moo  value  of 
0.95.  The  value  of  Ca  for  this  case  is  0.069,  which, 
as  expected,  falls  between  the  corresponding  values  of 
0.081  and  0.056  for  the  D  =  0.8  inch  and  D  =  1.6 
inch  separation  distances,  respectively. 

Figure  10  compares  the  relation  between  C &  and 
Mqo  ,  obtained  computationally,  to  the  corresponding 
relation  obtained  from  experimental  data.  Force  mea¬ 
surements  were  not  performed  in  the  experiments  de¬ 
scribed  in  Ref  (5).  The  pressure  data  measurements, 
however,  were  integrated  numerically,  by  Chapman 
and  Bonness,  to  produce  the  experimental  curve 
shown  in  Figure  10.  The  trends  between  both  curves 
are  similar,  however,  discrepancies  do  exist  between 
both  sets  of  data.  As  shown  earlier,  the  predicted 
pressure  coefficients  of  the  unfinned,  two  body  case 
were  shown  to  agree  well  with  experimental  data.  The 
relative  difference  between  the  experimental  interfer¬ 
ence  force  and  the  predicted  force  appears  large,  but 
the  absolute  difference  is  consistent  with  the  absolute 
differences  in  Cp.  Calculation  of  the  interference  force 
from  A Cp  for  this  body  of  revolution  results  in  small 
forces  and  a  magnification  of  the  relative  difference. 
This  is  true  for  both  the  experimental  and  computa¬ 
tional  approaches.  The  predicted  value  of  interference 
force  at  each  freestream  Mach  number  is  within  the 
measure  uncertainty  of  the  experimental  data.  The 
absolute  difference  in  the  forces  is  small,  and  for  a 
finned  body  with  larger  forces,  the  absolute  difference 
should  be  comparable  to  that  seen  here,  assuming  the 
flowfield  is  predicted  with  the  same ‘accuracy. 

In  comparing  the  computed  interference  force  to 
lift  forces  computed  from  experimental  data,  it  is 
found  that  the  interference  forces  computed  here  fall 
within  the  range  of  the  lift  experienced  by  a  single 
body  at  angles  of  attack  varying  between  4°  and  10°. 
This  estimate  should  only  be  viewed  as  an  approx¬ 
imate  one  in  view  of  the  discrepancies  observed  in 
Figure  10. 

The  three  store  mesh  described  above  was  used  to 
run  a  single  case  of  three  body  interference.  The  flow 
solution  was  obtained  for  the  parameters  Moo  =  0-95 
and  D  —  1.6  inches.  The  addition  of  a  third  body 
increases  the  flow  interference  level  in  comparison  to 
the  the  two  body  problem.  The  attractive  force  to¬ 
wards  the  axis  of  the  three  body  configuration  was, 
Ca  =  0.092.  This  represents  a  65%  increase  in  com¬ 
parison  to  the  two  body  problem.  Further  investiga¬ 
tions  of  three  body  interference  are  planned. 

CONCLUSIONS 

Flow  solutions  were  computed,  using  a  thin-layer, 
Navier-Stokes  solver,  about  interfering  bodies.  The 
computed  pressure  distributions  along  the  body  com- 
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pared  well  with  available  experimental  data  at  a  sep¬ 
aration  distance  of  1.6  inches. 

Three  regions  along  the  body  which  contribute  to 
repulsive  interference  forces  and  two  regions  which 
contribute  to  attractive  interference  forces  were  iden¬ 
tified.  The  combined  effects  of  these  regions  deter¬ 
mine  the  total  interference  force.  For  all  cases  inves¬ 
tigated  here,  that  force  was  an  attractive  force.  It  in¬ 
creases  with  free  stream  Mach  number  at  the  lower 
range  of  Mach  numbers  considered,  and  decreases 
with  free  stream  Mach  number  at  the  higher  range 
of  Mach  numbers  considered.  At  a  given  free  stream 
Mach  number  value,  the  interference  force  was  found 
to  increase  as  the  separation  distance  decreases. 

Flow  interference  measured  by  the  deviation  from 
symmetry  between  the  inboard  and  outboard  pres¬ 
sure  distributions  along  the  body  surface  show  that 
the  level  of  interference,  in  general,  increases  with 
free  stream  Mach  number  and  decreases  with  separa¬ 
tion  distance.  However,  at  the  upper  range  of  Mach 
numbers  considered,  the  increase  of  flow  interference 
with  free  stream  Mach  number  becomes  small  for  the 
large  separation  distance.  For  the  small  separation 
distance,  a  slight  reduction  in  interference  with  in¬ 
creased  Mach  number  is  observed  at  the  higher  range 
of  Mach  numbers  considered  here. 
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Figure  1(a).  Cross-Sectional  View  of  Two 
Store  Grid 


Figure  1(b).  Cross-Sectional  View  of  Three 
Store  Grid 


5 


Figure  1(c).  Lengthwise  Sectional  View  of  Two 
and  Three  Store  Grid 
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Figure  2(a).  Computational  versus  Experi¬ 
mental  Pressure  Distribution,  Two  Body  Case 


Figure  2(b).  Computational  versus  Experi¬ 
mental  Pressure  Distribution,  Two  Body  Case 


Figure  2(c).  Computational  versus  Experi¬ 
mental  Pressure  Distribution,  Two  Body  Case 


Figure  3(a).  Computational  versus  Experi¬ 
mental  Pressure  Distribution,  Two  Body  Case 
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Figure  3(b).  Computational  versus  Experi¬ 
mental  Pressure  Distribution,  Two  Body  Case 


Figure  3(c).  Computational  versus  Experi¬ 
mental  Pressure  Distribution,  Two  Body  Case 
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Figure  3(d).  Computational  versus  Experi¬ 
mental  Pressure  Distribution,  Two  Body  Case 


Figure  4(a).  Comparison  of  Pressure  Distribu 
tions  Showing  Interference  Force  Regions 
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Figure  4(b).  Aerodynamic  Interference  Com¬ 
parisons 


Figure  4(c).  Aerodynamic  Interference  Com 
pari  sons 


Figure  6(a).  Aerodynamic  Interference  Com-  Figure  5(b).  Aerodynamic  Interference  Com 
parisons  parisons 


Figure  6(b).  Variation  of  Cp  with  Increasing 
Separation  Distance,  D 


Figure  7.  Moo  Effect  on  Interference  Parame¬ 
ters,  D  —  1.6  Inches 


Figure  10.  Comparison  of  Computational 
and  Experimental  Interference  Force,  D  =  1.6 
Inches 


Figure  9.  Effects  of  MOT  and  D  on  Interference 
Force 
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